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Abstract 

We consider the resummation of soft-gluon emission in the transverse-mo- 
mentum distribution of vector mesons in hadronic collisions. We find that the 
resummed expression in the impact-parameter formulation has an expansion 
in a s with factorially growing terms with oscillating signs. These diverging 
terms arise from the small impact-parameter region of integration, and are of 
a subleading nature. We also obtain a closed expression for the next-to-leading 
logarithm resummation in g T -space, and we study its analytic structure. We 
find in this case that, although no factorially growing terms are present, there 
are geometrical singularities that severely restrict the range of applicability of 
the resummation formula. 
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1 Introduction 



In this paper we discuss soft-gluon resummation effects in the transverse momen- 
tum distributions of heavy systems produced in hadronic collisions. We will mainly 
focus on the Drell-Yan process, but several aspects of the problems we are considering 
here are of a general nature, and apply also to other cases, such as the production of 
Higgs bosons or of heavy-quark pairs. 

To simplify the notation, we consider the production of a virtual photon of off- 
shellness Q (this case can be trivially extended to W or Z production). We consider 
the kinematic region where its transverse momentum q T is much less than Q, but 
much larger than A QCD . This regime, characterized by two distinct scales, has been 
extensively studied in the past [[I], Large logarithms of the ratio Q/q T (the so- 
called Sudakov logarithms) arise to all orders of the perturbative expansion, with 
up to two powers of the logarithm for each power of a s . All-order resummation of 
the perturbative expansion is needed in order to make sensible predictions in this 
regime. In the small-g T limit, the resummed cross section has the structure suggested 
by Dokshitzer, Dyakonov and Troyan (DDT) in ref. [0: 

da Uq v-^ 2 d 



dQ 2 dq, 



ga ^ ] e g J dx A dx B 5{x A x B t) 

{f q /A(x A ,q T ) fg /B (x B ,q T )expS(Q,q T ) + (q ^ q)^ , (1.1) 



where a = Araa\ m l (9s) is the lowest-order cross section for the process qq — > ji + jj,' 
t = Q 2 /s, s is the hadronic centre-of-mass energy squared, and 



i(a s (/i 2 )) log % + B(a s (^ 2 )) 



(1.2) 



with 

A(a s ) = Axcis + A 2 c? s + . . . , B(a s ) = B ias + . . . . (1.3) 

In eq. ( |1 . 1|) , the leading logarithmic (LL) level of accuracy, which accounts for all 
terms of order a™l n+1 (where I = \og(Q 2 /q 2 )) in eq. (|1.2|) , is achieved by including 
only the A\ term of the expansion of A(a s ). The next-to-leading logarithmic (NLL) 
level, which amounts to the inclusion of all terms of order a™l n , is achieved if one 
includes also the A 2 and the B\ terms. Using q T instead of Q in the parton densities 
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also corresponds to a NLL correction.0 

We observe that no simple classification of the Sudakov logarithms can be given 
by looking at the expansion of the cross section ( |1 . 1|) in powers of a s and of I. This 
is because the exponent S has the formal expansion in a s 

oo oo 

S(Q, q T ) = <n+i<(Q 2 ) l n+l + Yl <na n s {Q 2 ) /" + .... (1.4) 

n=l n=l 

When only LL terms are included, the neglected NLL corrections correspond to a 
factor, which is a power series in terms of a s l, with the th order term equal to 1. 
Thus, as long as a s l < 1, the expansion is valid. On the other hand, by expanding 
the exponential of the first LL term, which is of order a s l 2 , we can get arbitrarily 
large powers of a s l 2 , and NLL corrections would appear to be of the order of a"/ 2n_1 . 
In this case, the expansion would seem to be valid only when a s l 2 < 1, a much more 
restricted range. 

As we will show later, eq. ( |1 . 1| ) has never appeared in the literature in a complete 
form at the NLL level. Indeed, it was shown in ref. @ that the resummation of 
soft-gluon effects is most naturally performed using the impact-parameter formalism. 
A general formula for the Drell-Yan cross section in the impact-parameter space has 
been obtained in ref. ||: 

d&dtf = £^J2<J d ^ *(*A* B - t) I d2b **** 

x {f q /A(x A ,ci/b) / 5/B (x B ,Ci/6)exp5(Q,6) + (q <-> q)} , (1.5) 



where 



S(Q,b) 



A(a s (// 2 ))log^f- + B(a s (^)) 



. 2 



,2 



;i.6) 



A(a s ) = A x a s + A 2 a z s + . . . , B(a s ) = B l0 c s + . . . , (1.7) 

and C\, C2 are arbitrary constants of order one.f] The coefficients A4 and Bi also carry a 
dependence upon these arbitrary constants (see |Appendix A|). There are well-known 



4 In order to go beyond the NLL terms, besides the inclusion of higher-order terms in the expansion 
of A(a s ) and B(a s ), it is also needed to perform a suitable change of scheme on the parton densities. 

5 In ref. ||, eq. ( |l.5[ ) is claimed to be valid also for q T < A QCD . This regime is not under 
investigation here. 



-3- 



problems that arise when applying eq. (1.5) to obtain phenomenological results. First 
of all, one can immediately recognize that, for any value of q T , eq. ( |1.5| ), as it stands, 
is undefined, because the b integration forces the scale at which a s and the parton 
densities are evaluated, to approach the non-perturbative region. This problem is 
usually avoided by replacing b in S(Q, b) and in the parton densities with M 



1 + {b/b lhQ f 



which approaches b for b <C &i im , and never exceeds bu m for large values of b. Choosing 
bn m of the order of 1 GeV -1 , one thus avoids the non-perturbative region. This pro- 
cedure introduces corrections of the order of powers of l/b\ im to eq. (|1.5[). Apart from 
the introduction of b*, eq. ( |1.5|) is usually also supplemented with a non-perturbative 
correction factorf] of the form exp(— gb 2 ). We stress that this factor is actually not 
needed in order for the integral in eq. ( |1.5| ) to converge (a proof is given in Appendix 



D 



In the present work, we will obtain an explicit expression for eq. ( |1.1|) , valid at 
the NLL level, and we will address several theoretical problems in the formulation of 
eqs. and QTJ). 



2 g T -space formulation of the resummed cross section 

In this section, we will obtain a g T -space expression for the cross section given 
in eq. ( |1.5| ), which is valid in the NLL approximation. We will analitycally perform 



the integration over 6, which appears in eq. ( |1.5|) . We will retain only the terms 
required for maintaining a NLL accuracy in the formal power expansion of the cross 
section. In other words, we will obtain an expansion for S(Q, q T ) that corresponds to 
retaining the first two sums of eq. (|1.4j) , and nothing more. We will first consider a 
simplified model, which corresponds roughly to the assumption of scale-independent 
parton densities and coupling constant. This model can be worked out quite easily, 
and it displays some features and problems that persist in the realistic case. 



6 In fact, it is argued in rcf. || that the non-perturbative region contributes a factor of the form 
exp[- S i(6) log(Q 2 )- 52 (&)]. 
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2.1 A simple model 

We thus assume that the parton densities and a s do not depend on the energy- 
scale. In this approximation, eq. ( |1.5| ) factorizes in a simple form: 

dQ^dq 2 = / dxAdxB fi/A(x A ) fq/ B (x B ) 5(x A x B -r), (2.1) 

where, following ref. 0, we have defined the form factor 

F(Q,q T ) = ^J d 2 be^ rb expS(Q,b). (2.2) 



After performing the angular integration in eq. ( |2.2j ) and changing the integration 
variable from b to b = bq T we get 

r\2 poo 

F(Q,q T ) = £L / dbb J (S) exp<S(Q,&/g T ), (2.3) 



where Jo is the th -order Bessel function of the first kind. Equation (|2.3| ) can be 
integrated by parts, exploiting the identity 

xJo ( x ) = A- (xJ^x)). (2.4) 
ax 

We obtain 

W,9t) = Q 2 TT / dbJ ^ b ) expS(Q,b/q T ). (2.5) 
In this way, we have an expression in the DDT form, analogous to eq. ( p..l[) . 
We assume that the Sudakov form factor at the LL level takes the form 

S(Q,b/q T ) = -aL\ (2.6) 

where 

L = \og^- (2.7) 

and 

a = ^-a s , (2.8) 



which is what is obtained if the running of a s is neglected in eq. ( |1.6| ). In this context 



we take the arbitrary constants ci, C2 to be equal to 1; different choices amount to a 
rescaling of the coupling constant. Defining 

Z = log^, ^ = log6 2 , (2.9) 
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we obtain 

F(Q,q, 



Q 2 -T~Z \ ex P(~ a ^ 

aq T 



db J\ib) exp (— 2all^ — a/~) 



(2.10) 



The exponential under the integral sign can now be expanded and integrated term by 
term. In fact, Ji(b) has an oscillating behaviour for large b, which gives rise to finite 
integrals when multiplied by any power of l b . Furthermore, all terms of the expansion 
will have at most as many powers of I as powers of a. The NLL level of accuracy is 
reached by neglecting the a/? term in the exponent, which is formally subleading. We 
thus get, at the NLL level, 



F(Q,q T ) = Q 2 



d 
dq\ 

d 

dq\ 



exp(— aV 



exp( 



-al 2 ) 2 



-Aal 



dbJi(b) b' 4al 

T(l - 2aiy 
T(l + 2a/) 



where the last integral was performed using the result 



dxx h Ji (x) 



F(l + h/2) 



-2 < h< 



1 



(2.11) 



f2.12) 



r(l -h/2)' ~ 2 

The form factor in eq. ( p. 11 ) is singular for al — > 1/2. Considered as a function of a at 
a given value of /, it is analytic in a circle of radius 1/ (2/) around a = 0. Thus, when 
expanded in powers of a, the coefficients of the expansion grow at most geometrically. 



On the other hand, the original expression eq. (|2.10|) is well defined for any value 
of a > and it is divergent for any negative value of a. Thus it cannot be analytic 
in the neighbourhood of the origin. This observation anticipates our conclusions: the 
original form factor has a divergent power series expansion in a, but, if we retain only 
the LL and NLL terms, we are left with a convergent expansion. Thus, the divergent 
behaviour must be ascribed to the subleading terms in the original expression. 

We now analyse more closely the origin of the non-analyticity of eq. (|2.10| ), now 
interpreted as defining a formal expansion in powers of a (it is easy to see that the 
coefficients of this expansion are finite order by order). First of all, we separate the 
large-6 and small-6 contributions to the integral. No singularities may arise from the 
large-6 region. This can be shown as follows. We rewrite the large-6 part of the 
integral in eq. fl2.10| ) as 



poo 2 roo 

/ dbJiip) exp (— 2all~ b — atf\ = Re / dbKi(-ib) exp (— 2all~ b — a/^) 

A. * " " n Ub 



[2.13) 
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where 60 is an arbitrary cut-off. We can now deform the integration contour counter- 
clockwise, letting 6 — > ib. Because of the asymptotic behaviour 

K&) * , (2.14) 

the integral is strongly convergent, with all its derivatives with respect to a, regardless 
of the sign of a. It thus defines an analytic function with infinite radius of convergence. 

We now study the small-6 contribution to the integral. Using the small-6 behaviour 
Ji(S) « 6/2, we get 

/•So l2 00 2n (r> \| 

/ dftSexp (-2o// g - all) = E^E^^^Tp + ^0)' ■ ( 2 - 15 ) 

J U 71=0 1=0 

We see that subleading terms in this expansion have factorially growing coefficients. 
For example, the terms corresponding to % = grow as (2n)\/n\. The factorial 
growth is accompanied by sign oscillations. This is why, when a > 0, this asymptotic 
expansion corresponds to a finite expression. 

In summary, we started with an expression that had an essential singularity at a = 
0, arising from the small-6 region of integration, but well defined for all positive values 
of a. By dropping all terms beyond the leading and next-to-leading, we obtained an 
expression that is analytic around the origin, but with poles at positive, finite values 
of a. We should now give our judgement, and decide which one of the two expressions 
is the most reliable. It is clear from the start that the factorially growing terms in 
our initial expression have no justification. It has been pointed out that, in fact, 
the expression for the Sudakov form factor should carry a theta function 

S(Q,b)^9(b-l/Q)S(Q,b), (2.16) 

since no large logarithms of b arise from the small-6 region in the exact calculation of 
the first-order correction to the Drell-Yan cross section. This theta function would 
suppress the integrand in the region < 6 < q T /Q. It would therefore amount to a 
power-suppressed correction that removes the small-6 singularities. In fact 

poo poo pQt/Q 

/ e- aL2e{b -^ /Q) J 1 {b)db= / e~ aL2 J x (b) db - 
Jo Jo Jo 

(2.17) 

If we expand the exponential in the first term, and integrate term by term, we find a 
power expansion in a whose coefficients are polynomials in log q T /Q. For the second 



-aL 2 



Ji(6) db . 
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term, the same procedure yields an expansion whose coefficients are suppressed by 
at least two powers of q T /Q, as can be seen by expanding the function J\ for small 
values of its argument. The factorial growth present in the second term compensates 
the one present in the first one. The fact remains, however, that if we drop all 
power-suppressed effects from our expression, the factorial growth remains. 

It should also be pointed out that the replacement given in eq. ( |2.16 ) is actually 
justified only for the first-order expansion of the Sudakov exponential. In fact, fac- 
torization of soft-gluon emission has a double origin. It comes from the dynamical 
factorization of the soft-gluon emission amplitudes, and from the kinematical factor- 
ization of the phase space for the emitted gluons. In the computation of the transverse 
momentum of DY pairs, the phase space factorizes naturally in the impact-parameter 
(i.e. in b) representation. However, neither the phase space nor the emission ampli- 
tudes factorize for the emission of hard gluons. To be more specific, while it seems 
reasonable to assume that gluon emission is cut off for gluons with transverse mo- 
mentum larger than Q (i.e. for b < l/Q), for two-gluon emission this constraint 
should become stronger, something of the form \q T1 \ + \q T2 \ < Q, and so on. It is 
therefore quite possible that the factorial growth coming from the small-6 region is 
only due to an improper treatment of the multigluon emission for large transverse 
momenta. Since the factorial growth is accompanied by sign alternation, it leads to a 
suppression of the small-6 region, thus giving a more stable result, which is however 
not fully justified. 

It is interesting to compare what we have found here with what was found in 
the case of resummation of threshold-enhanced Sudakov effects J7|. In that case, it 
had been found that factorial growth of the perturbative expansion was present in 
the x-space formulation of the resummation. Its origin is in the small-momentum 
region, and is caused by the fact that momentum conservation is violated in the x- 
space formalism. If one instead uses the Mellin transform (i.e. iV-space) formalism, 
which correctly implements momentum conservation, no factorial growth is present. 
Kinematic factorization of soft-gluon emission is present in this case in the iV-space 
formulation. It then turns out that no factorial growth is present when the resum- 
mation is performed in the most natural formalism. 

In the case we are considering here, instead, the factorial growth arises in the most 
natural formulation, which is the 6-space formulation. However, it does arise from the 
large-momentum region, which is unreliably treated by the resummation formalism. 
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It is generally believed that factorial growth in the perturbative expansion may 
arise from ultraviolet or infrared renormalons. In the present case, as in the case 
of ref. [J/j, the factorial growth does not arise from renormalons, since it is present 
also at fixed a s . In both cases, the factorial growth appears to be an artefact of the 
formulation. This is reassuring, since in most studies of power-suppressed effects it 
is assumed that no other sources of factorial growth exist besides renormalons and 
instantons. 



2.2 The realistic case 

The realistic case is more complicated, because the parton densities carry a de- 
pendence on the scale 1/b, which prevents from factorizing the cross section as in 
eq. ( |2.1| ). We therefore proceed as in ref. || and take the Mellin transform with 



respect to r (at fixed Q 2 ) of eq. ( |1.5|) : 

a N = I dTT N - 1 — -^-. 2 . (2.18) 



1 _ x Q 2 da 

<7 dQ 2 dq 2 



We thus find 



where 



cr N = A.J d 2 he^ rh C N { Cl /b) expS(Q,b) , (2.19) 



C , A r (/ i )=X^ e 9 / dXAX A 1 fq/A&A, AO / dx B X B 1 fq/ B (x B ,ll) . (2.20) 

q Jo Jo 
We now use the parton-density evolution equation to writeQ 

C N {c x /b) = C N (c 2 Q)ex P g N (Q,b) (2.21) 

with 

r4Q 2 du 2 

G N (Q,b) = -2 / J^ lN {a s {^ 2 )) . (2.22) 

The anomalous dimensions of the parton-density functions have a perturbative ex- 
pansion given by 

1n{u s ) = IxjfCLa + l2,Na 2 s + ■■■ ■ (2.23) 

We have therefore 



<TN = -^C N (c 2 Q)F N (Q, q T ), (2.24) 



7 In order to simplify the discussion, we only consider the non-singlet case here. 
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where 



F N (Q,q T ) = %- I d 2 be^ b exp<S(Q,6) expG N (Q,b). 



(2.25) 



The form factor F N (Q, q T ) can be computed in the same way as F(Q, q T ) in eq. (|2.2|). 
We thus obtain the analogue of eq. ( |2.5| ): 



F N (Q, q T ) = Q 2 



_d_ 

dq, 



dbJi(b) exj)S(Q,b/q T ) expQ N (Q,b/q T ). 



(2.26) 



T JO 



In this case, however, the integrand in eq. fl2.26j ) is not defined at large b. In fact, for 
large b the lower integration bound in eqs. ( |1.6j ) and ( |2.22| ) decreases and eventually, 
for some value of b, a s (fi 2 ) in the integrand is computed at the Landau pole. For the 
moment, we consider eq. (|2.26|) as a representation of its formal power expansion in 
a s (Q 2 ). The Landau pole contributes to spoil the convergence of this expansion. We 
will study its convergence properties in [Appendix C| . 



From eqs. (|1.6| ) and ( |1.7|) we can obtain an explicit expansion for S(Q,b), which 
turns out to have the form 



S(Q, b) = J2 < l fMs L), a s = a s (c 2 2 Q 2 ) , L = log 



i=0 



c 2 Q 2 b 2 
c\q 2 



(2.27) 



In the following we will consider the NLL expansion (in terms of L) , which corresponds 
to 



where 



S(Q,b) ~ S(a 8 ,L) = —f (a s L) + h{a s L) 



A 1 



fo(y) = -7T hv + ^og(l - boy)] 

°0 



(2.28) 



(2.29) 



h(y) 



b 2 



1 1 2/1 7 \ i b oy 

2 log (1 - b y) + 



+ 



log(l - b y) 
1 - b y 



A, 



h 2 



log(l - b y) + 



boy 



1 - b y 

and bo, bi are the first coefficients of the QCD (3 function: 

^^0^ = -b o a 2 s (l + b ias + ...) 



-i log(l - boy), 
b 



(2.30) 



_ 33-2n / 



1 153 - 19n f 
2tt 33 -2n f 



(2.31) 
(2.32) 
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Observe that fo(y) is of order y 2 , and fx(y) is of order y as y^O. We expect the 
remaining fi(y) to have a power expansion in y around y = 0. 

Equation (|2.27|) represents a logarithm expansion of the form factor S in terms 
of L. However, as previously observed, the large logarithm of our physical problem 
is I rather than L, which depends on the integration variable b. The leading and 
next-to-leading logarithm expansions of the form factor S in terms of I are 

1 



S LL (Q,b/q T ) 
<S NLL (Q,b/q T ) 



-fo(<x a l), 



— fo « s + fi(a, I) + — — jr- l~ h . 
a s d{a s I) 



(2.33) 



^2.34) 



On the other hand, Qn is a pure NLL term, since the scale dependence of parton 
densities gives rise to corrections, which have the form of a power series in a s l. Fur- 
thermore, from the definition given in eq. ( |2.22j ), we have 



Qn{Qi bj q T ) = Gn(Q, 1/<?t) + NNLL terms. 



(2.35) 



We first consider the LL approximation. In this case the b integration is trivial, 
because exp<S is independent of b at the LL level, as is clear from eq. (|2.33| ). The 
factor in eq. (|2.22|) being a pure NLL correction, we therefore find 



F N (Q,q T ) = Q 2 -^-expS ] 



dq\ 

which is, up to NLL terms, the result of ref. ||. 



-"LL ; 



We now turn to the NLL approximation. Using eq. ( 2.34j) we have 



(2.36) 



<>x]><S KLL (<9. b/q-y ) = (^| cxp3(<t. s ./). 



where 



h = 2 



df (a s I) 



1A\ boaj 



d(a s I) b 1 — b a s l 

Since at NLL is independent of b, we have 



d 



F N (Q,q T ) = Q 2 — 



expQ N (Q,l/q T ) exp§(a s ,/) ( — 

Ci 



(2.37) 



f2.38) 



db&j^b) 



Q 2 



dq\ 



expQ N (Q,l/q T ) exp§(a s ,/) 



2c 2 y T(l + h/2) 

c x ) r(i-/i/2) 



^2.39) 
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Equation (|2.39| ) represents our final result for the form factor Fn(Q,Qt) at NLL. It 
can be used to obtain an expression similar to eq. ( |1 . 1| ) for the cross section; in fact, 
using the procedure outlined in ref. ||, we may replace eq. ( [2.39| ) in eq. ( |2.19|) and 
obtain 



a N 



d 

dq, 

d 

dq. 



C N {c 2 Q) expQ N (Q, l/q T ) exp§(a s ,Z) 



2C2 
Cl 



r(l + h/2) 

r(i - h/2) 



C N (ciq T ) exp S(a 8 ,l) 

Ci 



2c 2 Vr(l + h/2) 



T(l - h/2) 



(2.40) 



where we have used eq. ( |2.21| ). The N dependence has been completely absorbed in 
the factor C7v(cig T ), so that the inverse Mellin transform can be performed immedi- 
ately and the cross section takes the form 



da 



d 



dQ 2 dq 



Xy 



T 



fq/A^A, ?t) fq/ B (x B , <? T ) exp§(a s , I) 



2C2 
Cl 



r(l + h/2) 

r(i-/i/2) 



+ (q++q)}, (2.41) 



where h is given in eq. Q2.38| ). Equation ( p.41| ) gives an explicit expression of the re- 
summed transverse-momentum distribution of DY pairs at NLL accuracy. This result 
bears some similarities with that of ref. 0, but is in fact different. The difference 
arises because we insist on keeping leading and next-to-leading terms according to 
the classification described in the Introduction, while in ref. || the emphasis is on 
finding a good appoximation to the 6-space result. A formula analogous to ours, but 
for the case of the energy-energy correlation in e + e~ annihilations, can instead be 
found in ref. |§. 

As in the case of the simple example treated in the previous subsection, our 
expression is analytic in a s for small a s , and displays a pole at h = —2, which 
corresponds to 

a, = - : 1 — . (2.42) 



l(Ax + b ) ' 

This can also be seen as a lower limit on the allowable values of q T : 



q T > Qexp 



1 



2a s (A 1 + b ) 



A 



QCD 



Q 



A 



QCD 



(2.43) 
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We thus find the disappointing result that our formula breaks down at a value of q T 
that is an increasing function of Q/A QCD . 

The full 6-space expression is instead well defined for all positive values of a s . 
However, it is not analytic around the origin. We will now prove that the small-6 
region contributions to the coefficients of its power expansion are at least as divergent 
as in the simple example examined in the previous subsection. We thus consider the 
integral 



f dbb exp{-^[a s b L + \og(l- a s b L)]\ (2.44) 
o la s »o J 



with L = \og(Q 2 b 2 / q T 2 ) , where, since we are considering the small-6 region, we have 
approximated Ji(6) ~ 6/2. We have fixed for simplicity 6o = 1. We will focus on 
the most subleading terms, that is to say, those carrying no powers of \og(q T /Q) in 
the expansion. These terms are obtained by simply taking q T = Q. Expanding the 
exponent in powers of a s , our integral becomes 

/ dbb exp 
Jo 




(2.45) 

where the highest power in h h comes from the exponentiation of the k = 2 term in the 
exponent, whose power of l~ b is twice the corresponding power of a s . The important 
thing to notice is that all the remaining coefficients r lt r 2 , etc., are positive. Thus 
each term of the sum is positive, and the value of each coefficient is bounded from 
below by the value of the first term db 6 (— /g) 2 - 7 = (2j)!/2, which shows the presence 
of factorial growth. 

In order to complete our task, we should show that no factorial growth arises from 
the large-6 region of integration 

/ d6Ji(S) exp j [a s b L + log(l - a s b L)\ 1 . (2.46) 

Unlike the case of the simple model, here we can no longer rely upon the analyticity 
in a s of this integral. In fact, the logarithm in the exponent encounters a Landau 
pole at a s b L = 1, and thus analyticity in a s cannot possibly be there. Thus, the 
coefficients of its perturbative expansion will grow faster than any power, but it is 
still possible to show that they grow less than factorially, by using the fact that the 



above integral can be related to an exponentially dumped one. The demonstration, 
similar to the analogous one in ref. J7[, is given in [Appendix C . 



3 Comparison of various approaches 



We will now compare the numerical differences of the various approaches for the 
computation of the Sudakov form factor. For the purpose of comparison, we will 



consider eq. ( |2.26| ) without the iV-dependent term (which is usually absorbed into a 
scale change in the parton densities) 



d 



F{Q, <? T ) = Q 2 I dbJx(b) ex P S(Q,V<fr 
aq T J 



(3.1) 



We will consider the following approaches: 



The full 6-space formula at NLL, that is eq. (|3.1|) with S(Q, b) given by 



and 



S(Q,b) = §(a s ,L* 
c 2 2 Q 2 b* 2 



-fo{a s L*) + h{a s V 



log 



42. ' 



(3.2) 



(3.3) 



The functions f and fi are given in eqs. ( |2.29| ) and ( |2.30| ). Here we choose &i im 
at the position that corresponds to the Landau pole. 



The NLL g T -space formula 



F(Q,q T )=Q' 



d 
dq\ 



exp S(a s , /) 



2c 2 \ h T{l + h/2) 



ci J T(l - h/2) 



(3.4) 



The expression of ref. 0. We wish to note at this point that the expression of 
ref. differs from our eq. (|3.4j ) by next-to- leading terms. In fact, the expression 
of ref. || is accurate at the leading-logarithmic level, with the further inclusion 
of subleading logarithmic terms down to the order a\l 2 ^ 2 in the expansion of 
the form factor. This is unlike our approach, in which terms down to the order 
a{P are kept in the Sudakov exponent. 
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Figure 1: Comparison of the numerical results of the various approaches for 
the calculation of the form factor. 



The results of the various approaches are summarized in fig. [I|. We see that the b- 
space, the NLL g T -space, and the EV approaches are quite similar for q T > 3 GeV. 
Below 3 GeV, the NLL g T -space formula begins to display a significant deviation, due 
to the fact that it is approaching the singularity h— > — 2. As an illustration, we also 
plot a modification of the NLL g T -space formula, obtained by letting 

5 NLL -> 6(b - 1/Q)5 NLL = 0(6 - g T /Q)<S NLL (3.5) 



before performing the b integration (the relevant formulae are given in |Appendix B|) . 



Modified in this way, we see that the NLL result is stabilized down to smaller values 
of q T . 

The result of ref. f| (the dashed curve labelled EV in the figure) is better behaved 
at small q T , and much closer to the full 6-space formula. In fact, the EV result differs 
from the 6-space formula by next-to-leading terms. It is paradoxically the lack of 
these terms (which would generate the singularity) that stabilizes the result. 
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4 Conclusions 

In summary, we have examined the impact-parameter formula for the transverse- 
momentum distribution in Drell-Yan-like processes. We have found that, when the 
resummed expression is expanded in powers of the strong coupling constant, and all 
terms suppressed by powers of q T /Q are neglected, the coefficients of the expansion 
exhibit factorial growth with oscillating signs. These diverging terms arise from the 
small-impact-parameter region of integration. 

The origin of these terms is analogous to what was found in the case of resum- 
mation of threshold-enhanced Sudakov effects 0. In that case, however, factorially 
growing terms arise in the x-space formulation of the resummation, from the small- 
momentum region of integration. No factorial growth is present in the iV-space formu- 
lation, which is most natural in this case, since it preserves momentum conservation. 

In the case we are considering here, instead, the factorial growth arises in the most 
natural formulation, which is the 6-space formulation. It does, however, arise from the 
large-momentum region, which is unreliably treated by the resummation formalism. 
Presumably, its origin is due to the fact that energy constraints are not correctly 
implemented by the resummation formula for multigluon emission of sizeable energy. 

In both cases, the factorial growth is only present in subleading (in the logarithmic 
sense) terms, which are not reliably given by the resummation formula. 

We have obtained an analytic expression for the resummed cross section in which 
only the LL and NLL terms in the Sudakov exponent are retained. Thus, factorially 
growing terms are not present in this formula. However, the formula displays geo- 
metric singularities at finite values of a s logQ/q T , which severely restrict its range 
of applicability. These geometric singularities arise from the small-6 region of inte- 
gration. In the original 6-space formula, the small-6 region gives rise to factorially 
growing terms, and thus to zero radius of convergence. This factorial growth is ac- 
companied by sign oscillation in this formulation, and thus by strong dumping of the 
small-6 region. Thus, the formula is well defined for positive values of the coupling 
constant. Since, however, the factorial growth is an artefact of the 6-space formula- 
tion, its superior stability (as far as the small-6 region is concerned) may hide some 
real problems. 

We believe that a more realistic treatment of the factorization properties of mul- 
tigluon emission at large transverse momenta could solve many of the problems dis- 
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cussed so far, and would thus be a valuable improvement of the computation of the 
Drell-Yan transverse-momentum spectra. 
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Appendix A Coefficients of the Sudakov exponent 

In the case of DY pair production, the coefficients in eq. (|1.7|) are given by [|l(| [11 1 



Ax — — (A.1) 

7T 



1 /67 vr 2 10 



57T c\e 



5i = ^log^- . (A.3) 

TT ZC 2 



The coefficient B 2 has also been computed: 



,67 TT 2 10 V cie^- 3 / 4 
2 T - T-n n f lQ g 



9 3 27 V 2c 2 
+ ^ 6 (log 2 - log 2 c 2 e 3 > 



9 7tt 2 2C.3 / 5 2tt 
~ 8 + "6" + T" + V 36 ~ ~W ' Uf 



(A.4) 



However, it only contributes at the NNLL level, together with unknown contributions 
from the A% term. 

Appendix B g T -space result with a b cut 

Our starting point is eq. ( |3.1|) , with 

S(Q, b/q T ) -> 0(b - p)S NLL (Q, b/q T ) , (B.l) 
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where p = q T /Q. We obtain 



d 



Q T2 / dhJ ^ h ) ex P d (b-p)S NLL (Q,b/q T ) 
aq T J 



Q 2 T^{ I db.h(b) l-ex V S Nhh (Q,b/q 



It K Jo 



+ / db Jx(b) exp«S NLL (Q,6/g 1 
Jo 

The second term in the curly bracket corresponds to eq. 
computed by expanding the Bessel function 

■M»> = E<W = 5-Ig + - 



k=l 



(B.2) 

The first term is easily 
(B.3) 



We obtain 

F(Q,q, 



> h n h+k+l 

-expS(a s ,l) ( — * ' 



Ci / /i + fc + 1 



(B.4) 



Notice that the term with k = 1 in the sum has a singularity for h = —2, which 
cancels exactly the analogous singularity arising in the second term from V(l + h/2). 
Similarly, the singularities for h = —4, —6 , . . . are cancelled by the higher-order terms 
in the sum. 



Appendix C Growth of the resummed expansion 

In this Appendix, we show that the large-6 region of the integral 

POO 

dbJ x {b)I{a s ,L) (C.l) 

Jb 

with 



I(a s , L) = exp 



1 



a s b 2 



[a s b Q L + log(l - a s b L)} 



(C.2) 



does not generate factorially growing coefficients in the perturbative expansion. 
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In the following, we will use the notation 

J2CjA-^y(-L) k <Y,G hk {-UsY{-Lt (C3) 

j,k j,k 

to mean that 

\C jt k\ < \G jtk \ (C.4) 

for all values of j and k. The minus sign in front of a s and L above is actually 
irrelevant, and it is there only to make the following discussion more transparent. We 
have 

T( n a s b L + log(l - a s b L) ~ ^^k-i^k 

/(«., L) = exp = exp £ . (C.5) 

s fc=2 U 

The coefficients of (— a; s ) (— L) k in the power series in the exponent are all positive, 
and the exponential function has a power expansion with positive coefficients. Thus, 
if we increase the coefficients of (— at s ) 1 (— L) k in the exponent, we will surely obtain 
an expression that is larger (in the sense of eq. (|C.4Q ) than the original one. We thus 
replace k — > k — 1 in the denominators of the coefficients, and obtain 

I(a s , L) < exp £ ° 1 h ( V ,\ ] = (1 " *sb Ly . (C.6) 

where y = L/b . We now continue with the trivial inequalities 

I{a s ,L) < (l-a s b L)y 

, h y{y-l)...(y-k + l) 



k\ 



k=0 

^ Y,^b Q L) k{ ^L < £(2a.&3)*(j/ + l) 2fc , (C.7) 

fc=0 fc=0 

where we have used the property 

k\ > k k 2~ k . (C.8) 
We thus have to examine the asymptotic behaviour of the integral 

POO 

+ 1) 2 *, (C.9) 

^6o 
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which in turn is that of an integral of the form 8 

db exp(-b) \og 2k b . (CIO) 



'60 

Using saddle-point methods, this is easily seen to grow like log 2fc 2k, which is stronger 
than geometric, but much less than factorial. 



Appendix D A few remarks on the non-perturbative part 

We consider here the 6-space resummation expression 
da do \ -» 



e 2 J dx A dx B 5(x A x B — r) J d 2 be l ^ T ' b 



dQ 2 dq% Q 2 ^ 

x {f q/A (x A ,ci/b*) fq/^x^Cx/b*) expS(Q,b*) + (q q)} , (D.l) 

(D.2) 



where 

v , i 

•A + &7C 

The new variable b* is approximately equal to b for small b, but never exceeds by m . 
We should have b\\ ni < 1/A QCD , so that both the Sudakov exponent (|1.6|) and the 
parton densities do not become undefined because of the Landau pole. 

The region of large b is in any case out of the control of perturbation theory, since 
it involves large values of the strong coupling constant. Therefore, one might want 
to keep b\ im substantially smaller than the value corresponding to the Landau pole, 
and parametrize non-perturbative effects by modifying the form factor in a way that 
leaves it unaffected at small b. A procedure that is often adopted is to multiply the 
integrand of eq. ( |2.3| ) by 

F NP = e- 9b \ (D.3) 

where the parameter g can be adjusted to the data. Notice that at large b, b* ap- 
proaches bn m , and therefore exp*S approaches a constant value while exp(ibq T ) oscil- 
lates. Nevertheless, the integral is convergent even in the limit g^O. In fact, if we 
add and subtract the contribution of exp S at b = &u m we obtain 

F(Q,q T ) = ^ J d 2 be^ S e- 9b2 [ex V S(Q,b*)-ex V S(Q,b hm )] 

8 Here we use the same method as was used in eq. ( |2.13[ ), and also use the fact that only the 
large-& region can lead to stronger-than-geometric growth. 
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+^e-^)ex P 5(g,fo lim ). (D.4) 

For g = the integral is now convergent at + oo for any choice of &i im , and for 
the remaining term we have 

lim ^! e -4/(4s) expS (Q, b hm ) = Q 2 n5 2 {q T ) expS(Q, b lim ) . (D.5) 

This fact may be useful to assess the relevance of the non-perturbative term compared 
to the effects of resummation. 
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